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Abstract 

It is shown that a single molecular magnet placed in a rapidly oscillating 
magnetic field displays the phenomenon of quenching of tunneling processes. 
The results open a way to manipulate the quantum states of molecular mag¬ 
nets by means of radiation in the terahertz range. Our analysis separates the 
time evolution into slow and fast components thereby obtaining an effective 
theory for the slow dynamics. This effective theory presents quenching of 
the tunnel effect. In particular, stands out its difference with the so-called 
coherent destruction of tunneling. We support our prediction with numerical 
evidence based on an exact solution of the Schrodinger’s equation. 
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1. Introduction 

Since a few decades molecular magnetic materials have arisen as a new 
test ground for several phenomena in quantum behavior of finite size mag¬ 
netic systems p El El a [5]. Molecular magnets have attracted attention 
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due to their potential in the implementation of several molecular spintronic 
devices [6l [TJ [8]. The long spin coherence times displayed turn them into 
promising candidates in the context of quantum computing, where the molecule 
spin is used to encode q-bitsfTUl El E]- 

The quantum mechanical degrees of freedom associated with molecular 
magnets can be manipulated and controlled with great accuracy by the ap¬ 
plication of external magnetic helds[2l El IH E]. It has been proposed [TO] 
that suitable manipulation of molecular magnets with time dependent mag¬ 
netic fields can be used to control the population of their quantum states 
thereby paving the road toward an implementation of a quantum computa¬ 
tion scheme known as the Grover’s algorithm. Within the same framework 
in this work we propose a way to control the quantum mechanical state 
of a molecular magnet by means of a rapidly varying magnetic field. We 
will show that radiation within the range of terahertz frequencies can be 
used to quench the quantum state of a molecular magnet thereby provid¬ 
ing a useful tool for potential applications. There is previous work relat¬ 
ing radiation in the terahertz range with magnetic properties of materials, 
for example terahertz radiation has been used to control the spin waves of 
antiferromagnets[IIl [I2] and in Ref. [13] hybrid magnetic structures have 
been used to generate radiation in the terahertz range. Our predictions al¬ 
low an extension of such control into the subject of molecular magnets. We 
study the dynamics of the magnetization of a molecular magnet exposed to 
circularly polarized terahertz radiation. In response to such perturbations 
the system can display a quenching of the tunneling rate between energeti¬ 
cally equivalent states giving rise to a trapping of the quantum mechanical 
state. This effect is analogous to an interesting effect in classical mechanics 
the trapping of a classical particle that can be achieved by introducing a 
rapidly oscillating potential V{x,t) = Vq{x) -\-Vi{x,t) [Tll[l5|. Under the ac¬ 
tion of such forces the slow dynamics of the particle is trapped by an effective 
potential Wfr = Vo{x) + /{2mLxP‘) where F is the force associated with the 

oscillating potential, F = —dxV\{x, t), uj its frequency and the bar represents 
an average over an entire cycle of the oscillating force. Based on those ideas 
it was proposed and demonstrated by KapitzajTS] that a pendulum with a 
rapidly vibrating point of suspension would be stabilized in the upward po¬ 
sition. Once stabilized the pendulum was shown to display small oscillations 
around its new equilibrium configuration. The main result of this paper is 
that a similar result holds for the quantum mechanical state associated with 
the spin of a molecular magnet. In this sense we can say that the prediction 
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corresponds to a Kapitza effect in Hilbert space. Instead of promoting tran¬ 
sitions between states, the high frequency radiation traps the state of the 
spin in a given configuration keeping it from describing tunnel transitions 
into other configurations. A quantitative statement of this effect is encoded 
in the tunneling time (the time that it takes for the system to tunnel from 
one minima to another) which is seen to diverge in certain circumstances. 
In response to the oscillatory disturbance the tunnel effect is suppressed and 
the states are frozen in a given configuration. A similar effect of suppression 
of tunneling has been reported in the literature concerning spins where the 
effect is attributed to interference of Berry’s phases of different paths asso¬ 
ciated with tunneling. In this context the effect has been dubbed quenching 
of the tunnel amplitude p6l l3]. 

2. Effective slow dynamics 

As is common in the theoretical studies of molecular magnets our study 
is based on the reduction of the electronic degrees of freedom, to some ef¬ 
fective low-energy Hamiltonian, based entirely on localized spin degrees of 
freedom of the magnetic ions within the molecule [IS]. In this context, the to¬ 
tal energy has two contributions H = 710 +Hi, one arising from the intrinsic 
anisotropy [2| , in the form: 

= -VSl + - S,2). (1) 

For T> > S this Hamiltonian represents a quantum spin with easy axis along 
the X direction and a hard axis along the .2 direction. The case £ = 0 
has been previously studied in Ref. ira This Hamiltonian can be used 
as a model for describing the magnetic degrees of freedom of several single 
molecule nanomagnets [2|. Among them the most widely used ones are the 
Mni 2 — ac moleculejS] (with S = 10, T> = 0.55 K and £ = 0.02 K), the 
molecular complex Feg molecular magnet[3| (with S = 10, T> = 0.29 K and 
£ = 0.05 K) and NR [23] (with S' = 4, I) = 0.75 K and vanishingly small 
£). The spin is perturbed by a circularly polarized time dependent external 
field [T2]: 

Hi = —h{cosu)t Sx + sin cut Sy) (2) 

Here cu is the frequency of the oscillation of the magnetic field in the x, y plane 
and h = gfisH with H being the amplitude of the oscillating magnetic field 
and g the gyromagnetic ratio (of order 1 in the examples given). Regarding 
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the geometry we can say that it is dominated by the anisotropy terms in the 
Hamiltonian. To apply our theory in the experimental setting the incident 
radiation must be polarized in the plane perpendicular to the hard axis. 
In a molecular magnet based crystal this can be achieved by selecting the 
crystal orientation with respect to the incident light. In the remaining parts 
of this article we will address the problem of how the spin responds to the 
perturbation in the limit of large to. In the special case of spin 1/2 all the 



Figure 1: Schematics of the proposed arrangement. A rapidly rotating magnetic field, 
h(t), in the plane induces a slow dynamics characterized by a stationary effective field, 
ffeff- The direction of the effective magnetic field is perpendicular to the plane. The slow 
dynamics is characterized by two degenerate classical ground states Cti and Cl 2 - In general, 
there are oscillations between those states mediated by quantum tunneling. However, for 
certain values of the parameters the quantum mechanical oscillations between those two 
states are quenched. 

anisotropic contributions reduce to the identity. The resulting behavior has 
the characteristic form of Rabi oscillations [ 21 ] . For higher spin the interaction 
between the magnetic moment and a radiation field is more complex and has 
been addressed in several references. For instance in Ref. [22| the effect of 
photon assisted tunneling events was reported in Feg samples irradiated with 
circularly polarized light. In Ref. [ 23 | microwave spectroscopy was used to 
reveal quantum superpositions of high spin states (S = 4) in NR. 

It is possible to derive a general treatment of a quantum system driven 
by rapidly varying potentials|22|. In general, the goal is to find an effective 
equation for the dynamics spanned by the Schrodinger equation: 

ihdt\'^) = H\^) (3) 

As in the classical case[25], the dynamics of the wave vector can be separated 
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into two components, a fast one that varies in a period of the potential and 
a slow one that evolves at a slower pace. The method consists in an exten¬ 
sion of the method of multiple scales from classical mechanics into quantum 
mechanics. We separate the slow and fast dynamics, by means of a unitary 
transformation, and proceed to write an effective theory that involves only 
the slow variables. The slow dynamics is affected by the rapid motion and 
is described by an effective Hamiltonian. We start from the time depen¬ 
dent Hamiltonian % and perform a time dependent unitary transformation 
exp(iW(t)). The basic idea is to absorb the time dependency of the Hamil¬ 
tonian into the operator W(t) and to obtain an effective time independent 
Hamiltonian given by: 

+ ih ( — j e-'^. (4) 


To find the specific representations of "Heff and W we write them as power 
series in l/o;: 


"Heff 




= n. 


( 0 ) 

eff 




1 


= —-I- - — 

nw {huy 


{huy 

+ 


^eff + 


( 5 ) 

( 6 ) 


Expanding Eq. 0 and equating both sides order by order we found recur¬ 
sive relations between and The expression for is found by 

enforcing that every term in the expansion of the effective Hamiltonian be 
time independent. A lengthy, but straightforward calculation leads to the 
following first contributions in the expansion for W: 


= —h {sin cotSx — cos cutsy) (7) 

and 


= -ih (cos ujt[Sx, 'Ho] + sin cot[Sy, Ho ]), 


while for the Hamiltonian we obtain: 


H 


( 1 ) _ 
eff ~ 



and 


H 


( 2 ) 

eff 
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{[[Sx,Ho],Sx] + [[Sy,Ho],Sy]) 


( 8 ) 

( 9 ) 

( 10 ) 
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Collecting the contributions up to second order we obtain an effective Hamil¬ 
tonian: 


"Heff = 'Hq — 


2huj 


S. 






i[[s,,no],s,] + [[Sy,Ho],Sy]), (11) 


In the limit of small VS/huj, i.e. for sufficiently rapid variations of the 
magnetic field, the effective Hamiltonian corresponds to the original station¬ 
ary contribution to the energy plus an effective field pointing in the direc¬ 
tion perpendicular to the plane of polarization of the magnetic field. For 
common molecular magnets such as the ones described above the condition 
VS/hjjj <C 1 sets the stage in a magnetic field oscillating in the terahertz 
rangejS]. The direction and magnitude of this effective field is independent 
of 'Ho • We can explain the origin of this contribution in terms of the following 
semi classical argument. The equation of motion of a spin in absence of any 
anisotropy is simply given by the Landau-Lifshitz equation, ^ = S x h{t), 
where h(t) correspond to the oscillating magnetic field. The motion can be 
understood as a succession of changes each of one is aprecession around the 
moving field. As shown in Fig. [^the complicated sequence of small preces¬ 
sions leads to a simple motion, a steady precession around the .2 axis. This 
can be associated with an effective field pointing along the same axis. As 


shown in Eq. (12) this field is generally present regardless of the base Hamil¬ 
tonian. This is not true for the following terms in the expansion. The second 
order term depends explicitly upon the stationary Hamiltonian. A funda¬ 
mental result of this paper follows after a direct calculation of this second 
order term, it can be easily verified that the contribution to order u)~^ merely 
shifts the values of the constants V and £ with a correction of order h?/ (/im)^. 
Hereafter we will denote those corrected values by T> and £ respectively. We 
conclude that the effective theory describes a spin with corrected anisotropy 
energies exposed to an effective magnetic field in the direction perpendicular 
to the plane of polarization of the time varying magnetic field. From this 
result we can infer several properties of the slow dynamics described by the 
quantum states. 


For spin 1/2 we can use our result in Eq. (12) in a direct way. Since the 
Hamiltonian "Ho reduces to the identity (a property of the spin 1/2 operators) 
the evaluation of the effective Hamiltonian is straightforward. We have a very 
simple result, namely: 
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Figure 2: Semiclassical origin of the effective magnetic field, (a) the spin precesses around 
an initial field hi. This precession is interrupted when the magnetic field changes into 
^ 2 - Under the new field the precession axes is changed and a new precessional motion is 
described. A second change into changes once again the precession axis. The net result 
of this sequence of changes is an effective precession around the ^ axis, (b) Exact solution 
of the Landau-Lifshitz equation for uj = bh. The exact results match the qualitative 
argument of panel (a). The results are consistent with a precession around an effective 
magnetic field pointing in the 2 ; axis. 


which provides a two level system with energy levels equal to A system 

prepared in a mixed state will oscillate [2T] between the eigenstates with a 
frequency determined by this splitting. For greater spins the anisotropy 
contribution plays an essential role in the dynamics as we will show in the 
next section. 

3. Dynamical suppression of quantum tunneling 

The slow dynamics reduces to the well-known problem of a spin under 
the action of an external magnetic field, along the hard axis, whose 

strength is given by = h?/{2huj). In this way we see that the resulting 

Hamiltonian is well studied in several contexts [2]. We can readily infer a 
number of properties regarding the behavior of the quantum moment by 
exploiting this analogy. We start by analyzing the classical limit. It is clear 
that the ground state is doubly degenerate. Using spherical coordinates, 
Cl = (sin 9 cos (/), sin 9 sin (/), cos 0), as shown in Fig. we find the two minima 
located at: 

COS 6*0 = - = - ^ 

4:Snw{V + £) 
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and 01^2 = 0, TT. The degeneracy of the classical ground states is lifted by 
quantum fluctuations associated with the tunnel effect between the two min¬ 
ima. The resulting energy splitting A between the two lowest lying states 
characterizes the tunneling time by the relation T = h/A. A detailed cal¬ 
culation of this splitting has been made and the results are highlighted in 
Fig. 4. A state initially prepared around one of the minima will oscillate, 
tunneling across the energy barrier, into the other minima in a characteristic 
time T. It is possible to And the tunneling t im e by direct diagonalization 
of the effective Hamiltonian. The tunnel splitting A oscillates as a function 
of the effective magnetic field. For certain specific values of the effective 
field a phenomenon known as the suppression or quenching of the tunnel 
effect can be observed. For those values of iifefr the tunnel splitting vanishes 
and the tunneling time diverges. The state is, therefore, trapped in a given 
state. The quenching of the tunneling processes, revealed by the reduction of 
the energy splitting between the two lowest lying energy states, has an ori¬ 
gin that can be traced back to the interference between different tunneling 
paths. Such behavior is better understood in terms of a semiclassical analysis 
of the tunneling process as is given by the instanton technique [26l [271 EH [29]. 
The origin of the oscillation of the tunneling gap is the interference between 
the Berry’s phases associated with complementary paths that accomplish 
the reversal. The result is[2S| A = AoCos(S'r2) where fl is the solid angle 
subtended by the complementary paths and Aq is a monotonous function 
of h. Whenever Sfl = {2n + l)n/2 for integer n the splitting vanishes. It 
is important to note that the dynamical supression of the Rabi oscillations 
due to the rapidly oscillating magnetic field arises due to the interference 
between tunneling amplitudes associated with different reversal paths and 
is, therefore, different in origin from the well-known coherent suppression of 
tunneling that has been predicted [30| in one dimensional systems. This ex¬ 
pression makes explicit the oscillations of the gap as function of the external 
driving field. However, the most important fact is that it emphasizes their 
origin in the Berry’s phase interference of complementary paths that achieve 
the tunnel reversal. For this reason, the effect is known as Berry’s phase 
interference quenching of the tunnel effect [26|. This semiclassical argument 
ought to be contrasted with the exact diagonalization of the Hamiltonian 
given in Eq. (12) for different values of the magnetic field intensity, h. After 
a direct numerical diagonalization, it is possible to calculate the energy gap 
between the two lowest lying states. This gap is shown, for <S' = l,5' = 3/2, 
S = 2 and S' = 5/2, in the continuous lines in the Fig. where the tunneling 
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times are plot as functions of h?/{2hjjj). In those figures it is evident the 
oscillatory behavior displayed by the gap and the quenching associated with 
its zeroes (indicated by arrows). We conclude that the effective slow dynam¬ 
ics of a spin under a circularly polarized magnetic field will be described by 
oscillations between two states |r2i) and 1 ^ 2 ), where r 2 i _2 correspond to the 
unit vector with polar angles (6q, 01,2)- Here, the state |r2) corresponds to the 
spin coherent state oriented along the direction of fl[T6j. Those oscillations 
correspond to tunneling events. Furthermore, it can be concluded that those 
oscillations are quenched for certain values of the amplitude of the oscillating 
magnetic field. In the next section this prediction will be contrasted with 
exact numerical results. 

4. Comparison with exact results 

We now contrast the predictions made so far with exact numerical results. 
To that end we solve numerically the time dependent Schrodinger equation 
Eq. If. To characterize the slow dynamics generated by the Schrodinger 
equation we have prepared the state initially in the ket |r2i) and let it evolve. 
We have computed the projection of the state vector |d^(t)} on the kets |fli) 
and \fl 2 ) with typical behavior shown in Fig. where we plot |(ni|\I/)p 
and |(n 2 |\k)p. We see that the projection on the initial state l^i) starts 
at its maximum value, one, and then it is reduced to zero. Meanwhile the 
projection on the second state describes the opposite behavior. We remark 
that the oscillation in the amplitude is accompanied by a rapidly oscillatory 
component that corresponds to the fast dynamics. This behavior, reminiscent 
to that of a two level system, corresponds to the quantum tunneling between 
the states described by the states |r2i) and 1^22) • The characteristic time T 
that takes a transition from one state to the other correspond to a direct 
quantification of the tunneling time. This time can be inferred directly from 
the numerical solution. For different values of the field we calculate such 
tunneling time. The results, for 5 = 1, S' = 3/2, S = 2 and S = 5/2, are 
displayed in Fig. 

In the different calculations we have chosen units in which I) = 1 and se¬ 
lected a value for S = 0.1T>. For the numerical calculations we have selected 
the driving frequency as a; = bOV/h. The results are in evident agreement 
with those obtained from the effective theory. This provides strong evidence 
in favor of the validity of the effective Hamiltonian (Eq. PH) to describe 
an important aspect of the complicated dynamics displayed by the system. 
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Figure 3: Characteristic behavior of the projections |(f2i|5')p and |(f22|^')P as functions 
of time. The initial state is |ffi) and in the slow dynamics regime the system manifest 
well defined tunneling oscillations between the two classical ground states. The tunneling 
time, defined as the characteristic time that takes the system from one configuration to 
the other, is denoted by T. The slow dynamics is dressed by very fast oscillations. 

These exact results show that the tunnel effect is effectively quenched by the 
time dependent magnetic field for certain values of the field intensity and 
frequency. This quenching is made manifest by a divergence of the tunnel¬ 
ing time associated with quantum transitions between the two degenerate 
classical ground states. 

5. Discussion 

In this paper we have presented a detailed study of the dynamics described 
by a quantum spin when exposed to a rapidly varying magnetic field. The 
problem is analogous to the pendulum, proposed by Kapitza, with a rapidly 
oscillating suspension point. Just like in the Kapitza’s pendulum the motion 
of the quantum states is decomposed into two components. We have a rapid 
behavior whose characteristic time is given by the external potential and a 
slow contribution. The slow dynamics is affected by the rapid motion and 
is described by an effective Hamiltonian. By performing this separation in 
a systematic fashion we have proved that the slow dynamics is accounted 
for by an effective Hamiltonian that has a Zeeman-like contribution in the 
direction perpendicular to the plane of polarization of the magnetic field. 
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Figure 4: Tunneling time is shown for spin S = 1, S = 3/2, S = 2, and S = 5/2 as 
a function of the intensity of the magnetic field hg// = In? The tunneling time, 

obtained directly from the numerical solutions of Eq. is represented by the dots. It is 
compared with the predictions from the effective theory given by Eq. (12). The effect of 
quenching of the tunnel effect implies a divergence of the tunneling time as indicated by 
the arrows. 


The effective Hamiltonian displays the interesting effect of geometric-phase 
induced quenching of the tunnel effect that can be used to freeze the quantum 
states in well defined configurations. Just like the Kapitza’s pendulum is 
trapped in the vertical position, unstable when the driving force is absent, 
the spin of the molecular magnet is trapped in a given quantum state. In this 
sense we talk about a Kapitza effect in Hilbert space. This effect provides a 
tool to control the quantum states by means of high frequency fields. For a 
typical molecular nanomagnet our analysis sets the desired frequencies in the 
terahertz range. In addition we have presented numerical evidence, based on 
a direct solution of the Schrodinger equation, that confirms the predictions 
of the effective description. 
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